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ABSTRACT 

We present age estimates for the newly discovered very r-process enhanced metal-poor star 
HE 1523-0901 ([Fe/H] = -2.95) based on the radioactive decay of Th and U. The bright {V = 11.1) 
giant was found amongst a sample of bright metal-poor stars selected from the Hamburg/ESO sur- 
vey. From an abundance analysis of a high-resolution [R — 75,000) VLT/UVES spectrum we find 
HE 1523—0901 to be strongly overabundant in r-process elements ([r/Fe] = 1.8). The abundances of 
heavy neutron-capture elements [Z > 56) measured in HE 1523—0901 match the scaled solar r-process 
pattern extremely well. We detect the strongest optical U line at 3859.57 A. For the first time, we are 
able to employ several different chronometers, such as the U/Th, U/Ir, Th/Eu and Th/Os ratios to 
measure the age of a star. The weighted average age of HE 1523—0901 is 13.2 Gyr. Several sources of 
uncertainties are assessed in detail. 

Subject headings: Galaxy: halo — stars: abundances — stars: individual (HE 1523—0901) — early 
universe — nuclear reactions, nucleosynthesis, abundances 
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1. INTRODUCTION 

The Galactic chemical history can be reconstructed by 
studying the most metal-poor stars. Of particular im- 
portance are direct age measurements of the oldest halo 
objects. Such ages can be inferred by comparing the 
observed abundance ratios of the radioactive elements 
232Th (half-life of 14 Gyr) and ^^^\J (4.5 Gyr) with the- 
oretical pred i ctions of their initial production. In 2001, 
ICavrel et al.l ()200ir ) announced the first detection of U 
in an r-process enhanced metal-poor star, CS 31082- 
001 ([Fe/H] = -2.9). Together with a Th measure- 
ment, this led to the novel use of the U/Th ratio as a 
chronometer. Th is technique yielded an age of 14 ±3 Gyr 
(jHill et al.ll2002t l. in good agreement with the recent re- 
sult of 13.7 ± 0.2 Gyr for the age of the Universe , as de- 
termined from WMAP data (jSpergel et al.ll2006f) . 

The similarities of the neighboring nuclei Th (Z=90) 

Based on observations made with ESO telescopes at the La 
Silla Paranal Observatory under programme ZD's 275.D-5028 and 
077.D-0453. 
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and U (Z=92) causes this cosmic clock to be 
much less prone to systematic uncertainties than the 
more widely used Th/Eu (Z— 63 for Eu) ch ronometer 
(|Gorielv fc Clerbauxlll999l : ISchatz et al.ll2002[ ). Hence, it 
provides more reliable age estimates. It is thus important 
to attempt a U measurement in further metal-poor stars. 
The strongest U line at optical wavelengths is very weak, 
requiring very high quality spectroscopic data. Addition- 
ally, it is located in the wing of a strong Fe line and is 
blended with a CN feature. Hence, low C and N abun- 
dances are vital for an accurate U abundance determi- 
nation. It is thus not possible to measur e U in C-rich 
r-prq cess rich stars such as CS 22892-052 ([Sneden et al.l 
Il996f ). for which only the Th/Eu chronometer can be 
employed. These observational difficulties explain why, 
currently, CS 31082-001 is the only known strongly r- 
process enhanced met al-poor sta.r with a reliably mea - 
sured U abundance fCavrel et al.l 1200 ll : iHiTl et al.ll2002f ). 
However, one further object is known to exhibit U (V. 
Hill et al. 2007, in preparation), while for a nother one, at 
least a tentative detection has been made ()Cowan et al.l 
[200I . 

Here we report the discovery of a strongly r-process 
enhanced star with a U measurement. This detec- 
tion, in combination with other neutron-capture abun- 
dances, makes possible stellar nucleo-chronometry. The 
giant HE 1523-0901 {V = 11.1) was found in a sam- 
ple of bright metal-poor stars (Frebcl et al. 2006) from 
the Hamburg/ESO Survey (HES). Based on an initial 
metallicity estimate of [Fe/H] = —2.7, obtained from a 
medium-resolution {FWHM ~ 2A) spectrum, the star 
was observed with high spectral resolution (i? ^ 45, 000) 
with MIKE at the Clay Magellan telescope. The spec- 
trum revealed a strong enhancement of heavy neutron- 
capture elements, and a tentative detection of U was 
made. 
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To confirm the tentative U detection in HE 1523—0901, 
new observations were obtained wi th the Ultrayiolet- 
Visual Echelle Spectrograph (UVES; iDekker et aLll2000D 
at the ESO Very Large Telescope. The star was observed 
in service mode in 2005 August and 2006 April. We 
made use of the image sheer No. 2 and a 0'.'45 slit width 
to achieve a very high resolving power of i? ~ 75, 000. 
We employed the BLUE 346 nm setting covering 3050- 
3874 A and the BLUE 437 nm setting covering 3758- 
4990 A. The two settings conveniently overlap at 3860 A 
where the U line is located. The total exposure time was 
7.5 h, of which 5 h were spent with the BLUE 346 nm and 
2.5 h with the BLUE 437 nm setting. 

We use the pipeline reduced spectra provided by the 
ESO Data Management and Operations Division. The 
frames, shifted to the stellar rest frame, are co-added. 
We estimate a signal-to-noise ratio of S/N ^ 350 per 
12.4mA pixel at - 3900 A, the region where the two 
settings overlap. 

3. ABUNDANCE ANALYSIS 

For our ID LTE abundance analysis of the VLT/UVES 
spectrum we use the latest version of the MARCS code 
(B. Gustafsson et al. 20 07, in preparat i on). Solar 
abundances are taken from lAsplund et al.l ()2005[ ) . For 
the choice of atomic absorption lines of the lighter el- 
ements we use a line list based on the com pilations of 
iBarklem et"all (|2005f ) and lAoki et al.l ()2002[ ) as well as 
our own collection retrieved from the VALD database 
(iKupka et al.lll999l ). The fine list provided in iHih et all 
(|2002f ) was updated with the latest atomic data used for 
the measurements of the neutron-capture elements. The 
molecular line data for CH (B. Plez, private communica- 
tion) is based on gf- values and line positions from LIE- 
BASE (Luque & Croslev"1999'), excitation energies are 
taken from Jorgcnscn et al. (199_6.), and isotopic shifts 
are computed by Plez. The CN lin e data (B. Plez, p ri- 
vate commu nication) is de s cribed in I Hill et al.l (|2002f ). 

Using the lAlonso et al.l (|1999D calibration, we deter- 
mine an effective temperature of T^s = 4630 ± 40K 
(random error) from dereddened BVRI CCD pho- 
tometrv (iBeers et all |2006D and JHK 2MASS data 
(jSkrutskie et al.l [20061 ). The microturbulence was ob- 
tained by demanding no trend of elemental abundances 
with equivalent width. Using Fei lines, the result is 
Vjiiicr = 2.6 ± 0.3 km s^^. From the LTE Fei-Feii ioniza- 
tion equilibrium we derive a surface gravity of logg = 
1.0±0.3. The resulting metallicity is [Fe/H] = -2.95±0.2 
for HE 1523-0901. 

Equivalent widths are obtained by fitting Gaussian 
profiles to the chosen atomic lines. For blended lines 
and molecular features we use the spectrum synthesis 
approach, in which the abundance is obtained by match- 
ing the observed to a synthetic spectrum of known abun- 
dance. This technique is extensively used for the Thii 
line region around 4019 A, the Uii line region around 
3860 A and other important line regions (e.g., Eu, Ir, Os). 
Abundance uncertainties arising from this method are 
usually driven by difficulties with the continuum place- 
ments and vary from 0.05 — 0.15 dex. 

4. NEUTRON-CAPTURE ELEMENTS 

The spectrum of HE 1523—0901 shows numerous 
strong lines of ^ 25 neutron-capture elements associ- 



ated with the r-process. In Figure [TJ the observed abun- 
dance pattern of the star is shown. Error bars repre- 
sent the standard error of the mean abundance of sev- 
eral lines for each element. As can be seen, it closely fol- 
lows the sca l ed ([r /Fe] = 1.8) solar r-process pattern of 
I Burr is et al.l (|2000( ). For neutron-capture elements with 
56 < Z < 77, the agreement is excellent (standard error 
of the mean is 0.02 - see bottom panel of figure). 

Several Th ii lines are detected. However, most of them 
are severely blended with lines from other elements. We 
regard the line at 4019.13 A the most reliable one be- 
cause the contaminating blends are best known from the 
literature. First, a strong ^'^CH feature is located blue- 
ward (4019.01 A) of the Th line. Despite the star having 
i^C/i^C of ~ 3-4, (based on features at - 4020, - 4220 
and 4310 A), the total carbon abundance is subsolar, 
[C/Fe] — —0.3, and thus not a major contaminator. Sec- 
ond, the blending with the Ceii line at 4019.06 A can 
easily be accounted for by fitting the blue wing of the ob- 
served feature. Fortunately, the red wing is dominated 
by the Th line and allows a well constrained fit. Any 
additional contamination of the whole region with nu- 
merous Co I lines is not significant. 

Other strong Thii lines are located in bluer regions of 
the spectrum with lower S/N. Also, the line fits suffer 
from unidentified features and blends that are difficult 
to account for. Due to these problems, we adopt the 
abundance of loge(Th) = -1.20 from the Thii4019A 
line. 

In our high-resolution, high S/N spectrum of 
HE 1523-0901, we detect the Uii line at 3859.6 A. The 
spectral region around this line is shown in Figure [2] We 
repeat here that low C and N abundances are crucial for a 
U detection because it is significantly blended with a CN 
feature. Both the U line and the CN feature, in turn, are 
located in the wing of a very strong, saturated Fei line. 
This complicates the fitting procedure. To obtain a best 
fit of this region^^, we increase the VALD log 5/ value of 
the Fei line at 3859.911 by 0.24 dex. The C ([C/Fe] = 
—0.15) and N ([N/Fe] = 0.6) abundances employed in 
the U synthesis are well constrained by a good fit to the 
violet CN feature at 3883 A. The U abundance is then ob- 
tained by simultaneously fitting the Fe, CN, and U lines. 
We derive an abundance of loge(U) — —2.06. The final 
U/Th ratio for HE 1523-0901 is loge(U/Th) = -0.86. 
Further abundance ratios of radioactive to naturally oc- 
curring r-process elements are listed in Table [T] 

The abundances of Eu, Os, and Ir are loge(Eu) = 
-0.62 loge(Os) = 0.18, and loge(Ir) = 0.24. We also 
attempted to detect Pb, the decay product of Th and 
U. However, it could not be detected in the current spec- 
trum of HE 1523-0901 (the S/N in this region is - 150). 
As shown in Figure[l] the upper limit of loge(Pb) < —0.2 
is below the expected abundance of the scaled r-process 
pattern. A full discussion of the complete abundance 
analysis will be given elsewhere (A. Frebel et al. 2007, 
in preparation). 

To test our derived abundances, we measured the Thii 

We note that the two features in the red wing of the Pel at 
3859.6 A line (see Piguro [Jjl were previously unidentified (Hill et 
al. 2002). Based on the new Sm atomic data l|Lawler et al.ll2"006l) 
we find the bluer feature (3860.28 A) to be a Sm line. The other 
feature remains unidentified. 
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A 4019 and the U feat ures in the spect rum of CS 31082- 
001 that was used bv iHill et all ()2002f ). Figure [2] shows 
the U region for CS 31082-001 (crosses). Despite differ- 
ences in the employed model atmospheres, we obtain a 
loge(U/Th) ratio of -0.93 for CS 31082-001. This is in 
very good agreement with the pubhshed value of —0.89, 
as derived from these two lines. 

We estimate a fitting uncertainty of 0.05 dex for the Eu, 
Os, Ir, and Th abundances. The U abundance is driven 
by the fit of the Fe line close to the U line. Changing 
the C by -fO.ldex results in only a —0.02 dex different 
U abundance. Changing the Fe abundance by -1-0.1 dex 
changes the U abundance by —0.12. We adopt a 0.12 dex 
uncertainty for U. 

5. NUCLEO-CHRONOMETRY 

There are three types of chronometers that involve the 
abundanc es of Th, U and na turally occurring r-process 
elements (jCavrel et all |200H) . The subscript "initial" 
refers to the initial production ratio (PR), while the sub- 
script "now" refers to the observed value. 

1. At = 46.7[log(Th/r)i„itiai - log e(Th/r)now] 

2. At = 14.8[log(U/r)i„itiai -loge(U/r)now] 

3. At = 21.8[log(U/Th)initiai - loge(U/Th)„ow] 

Using several different chronometers and PRs, we de- 
rive a set of ages for HE 1523—0901. The results are 
given in Table [1] Where available, we list several PRs 
for each chronometer to illustrate the available range and 
the subsequent spread in the derived ages. We take the 
weighted average of all the individual ages to derive a 
final age of 13.2 Gyr for HE 1523-0910. 

Forming an average based on weights obtained from 
the uncorrelated observational uncertainties is an arbi- 
trary choice which only minimizes the observational (sta- 
tistical) uncertainties, but not necessarily the systematic 
ones. Using different weights, for example by omitting 
the Th/r ratios, would lead to slightly larger observa- 
tional but smaller systematic uncertainties. A weighted 
observational uncertainty in the abundance ratios arising 
from the fitting procedure results in an 0.7 Gyr weighted 
uncertainty for the final age. This value is driven by the 
uncertainty of the uranium abundance measurement. 

We also investigate the influence by variations of model 
atmosphere parameters (Toff, log 5, Wmicr) on the the stel- 
lar age. Combining (square-root of quadratic sum) these 
three age uncertainties yields a 1.5 Gyr weighted uncer- 
tainty in the final age. Any correlations of the different 
chronometers are thus automatically taken into account. 
To obtain an age uncertainty arising from the uncertain- 
ties in the PRs we calculate at = J^ii'^i ^ '^^i)/ 
(with Wi = 1 /erf ; ati is the age uncertainty from the dif- 
ferent PRs and ai the one from the observational uncer- 
tainty) as an upper bound, assuming the worst possible 
correlation(s) of the uncertainties in the PRs. We thus 
derive a 2.7 Gyr weighted uncertainty in the final age. 
For the calculation of the PR uncertainties we followed 
Schatz et al. (2002) who list overall systematic uncer- 
tainties for all three types of chronometers. In Table [U 
we list the five age uncertainties for all chronometers. 

Due to the much shorter half-life of U, uncertainties 
in ages derived from chronometers U/r are significantly 



smaller than for those from Th/r. Reducing the num- 
ber of available chronometers by the Th/r ones yields a 
weighted average of 13.4 Gyr. The observational uncer- 
tainty becomes 0.8 Gyr and the combined model atmo- 
sphere uncertainty is 0.9 Gyr, while the PR uncertainty 
is much reduced to 1.8 Gyr. This illustrates the superior- 
ity of the U/r ages in terms of systematic uncertainties. 
Also, ~ 2 Gyr is roughly the age spread caused by the 
different PRs for U/Th. 

Another advantage of the special case of U/Th is that 
systematic observational uncertainties associated with 
model atmosphere parameters are less severe than for 
other ratios because the neighboring U and Th nuclides 
have very similar atomic parameters. It follows that 
the U/Th PR should be less sensitive to theoretical r- 
process model uncertainties compared with ratios of ele- 
ments with larger mass separation such as the case of T h 
to Eu (|Gorielv fc Clerbauxl [1991 iWanaio et ai][200l 
Howev er, in their extensive parameter studv lSchatz et al.l 
(|20f)l find very similar overall systematic uncertainties 
for the three chronometer types. 

In order to derive reliable ages for individual metal- 
poor stars it is crucial to have available stars with a very 
precise U measurements together with the most accurate 
theoretical PRs. Only then can nucleo-chronometry be 
refined to deliver the best possible Galactic stellar ages. 

6. DISCUSSION 

For the first time it is possible to make use of all 
three types of chronometers in one star: Th/r, U/r and 
U/Th. Previously, either the chronometer Th/r was em- 
ployed (e.g. CS 22892-052: ISneden et al.lll996D or only 
the ones based on the availability of U (CS 31082-001; 
ICavrel et al.ll200ll ). Also, the U/r chronometers, such as 
U/Os or U/Eu, for the first time deliver reasonable ages 
because, unlike CS 31082-001 (see Schatz ct al. 2002), 
HE 1523—0901 appears to follow the solar r-process pat- 
tern from Ba all the way to U. Given the observational 
and theoretical uncertainties, the consistency of the in- 
dividual ages derived for HE 1523—0901 is remarkable. 

Since we also measure the U/Th ratio in CS 31082- 
001 it is possible to determine a relative age of the two 
stars. We find HE 1523-0901 to be 1.5 Gyr younger 
than CS 31082-001, which is independent of the em- 
ployed production ratio. This age difference is based 
on only a 0.07 dex difference in the observed U/Th ra- 
tios. Given that the observational uncertainties exceed 
that ratio difference as well as the uncertainties associ- 
ated with initial production ratios, the present ages of 
the two stars suggest that they formed at roughly the 
same time. This is also refiected in their almost identi- 
cal metallicity. We note that HE 1523-0901 is - 200 K 
cooler than CS 31082-001, which makes the lines in our 
star more prominent (see Figure [2]). 

Stellar age measurements such as these provide a lower 
limit to the age of the Galaxy and hence, the Universe. 
Despite their large uncertainties the age limits provided 
by HE 1523-0901 and CS 31082-001 are in good agree- 
ment with the WMAP result of 13.7 Gyr for the age of 
the Universe. 

We note here that HE 1523—0901 is the first "uranium 
star" discovered in the HES. Only a very few r-process 
enhanced stars are suitable for a detection of U because 
the objects need to be bright, sufficiently cool, strongly 
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overabundant in heavy neutron-capture elements, and 
have low C and N abundances. It is of great impor- 
tance, however, to find further "uranium stars" . This 
new group of objects will provide crucial observational 
constraints to the study of the r-process and its possible 
production site(s). They may also be used to empirically 
constrain the U/Th production ratio, which, in turn, may 
provide the much-needed feedback for theoretical works 
to improve the initial production ratios. 
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TABLE 1 

Ages derived from different abundance ratios 



X/Y 


log(PR)'' 


Ref. 


loge(X/Y)obs 


Age (Gyr) 


Uncertainties (Gyr) 


Th/Eu 


—0.377 


1 


-0.58 


9.5 


3.3/3.4/0.6/0.6/5.6 


Th/Eu 


-0.33 


2 


—0.58 


11.7 


3.3/3.3/0.5/0.5/5.6 


Th/Eu 


-0.295 


3 


-0.58 


13.3 


3.3/3.0/0.2/0.2/5.6 


Th/Os 


-1.15 


2 


-1.38 


10.7 


3.3/2.8/5.6/0.0/5.6 


Th/Ir 


-1.18 


2 


-1.44 


12.1 


3.3/1.9/2.8/1.4/5.6 


Th/Ir 


-1.058 


1 


-1.44 


17.8 


3.3/2.0/2.9/1.5/5.6 


U/Eu 


-0.55 


2 


-1.44 


13.2 


1.9/0.6/0.4/0.2/1.6 


U/Os 


-1.37 


2 


-2.24 


12.9 


1.9/0.6/1.2/0.3/1.6 


U/Ir 


-1.40 


2 


-2.30 


13.3 


1.9/0.3/0.3/0.7/1.6 


U/Ir 


-1.298 


3 


-2.30 


14.8 


1.9/0.3/0.3/0.8/1.6 


U/Th 


-0.301 


4 


-0.86 


12.2 


2.8/0.4/0.9/0.4/2.2 


U/Th 


-0.29 


5 


-0.86 


12.4 


2.8/0.4/0.9/0.4/2.2 


U/Th 


-0.256 


3 


-0.86 


13.1 


2.8/0.5/1.0/0.5/2.2 


U/Th 


-0.243 


6 
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^ PR: initial production ratio^ Age uncertainties arising from uncertainties in : 
observ. measurements, T^ff, log 5', t'micn 

PR Ref s. — 1: ISneden et al.l 112003) . 
2: Schatz et al. (2002), 3: Cowan ct al. ( 200l) . 4: IGorielv Arnouldl Il200ll )~5: 
IWanaio et al.. L200%] . 6: ,Dauphas. (,2005i) 
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Fig. 1. — Neutron-capture element abundances of HE 1523—0901 in comparison with those from the solar r-process (jBurris etaDllQOO) 
scaled to match the observed elements with 56 < Z < 69 [top panel). The bottom panel shows the residuals from the abundances of 
HE 1523—0901 minus the solar r-process values. 
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Fig. 2. — Spectral region around the U II line in HE 1523—0901 {filled dots) and CS 31082-001 (crosses; right panel only). Overplotted 
are synthetic spectra with different U abundances of loge(U) = none, —1.96, —2.06, and —2.16 (HE 1523—0901) and loge(U) = none, 
—2.05, —2.15, and —2.25 (CS 31082-001). The dotted line in the left panel corresponds to a scaled solar r-process U abundance present in 
the star if no U had decayed. Positions of other features are indicated. 



